Josephson tunneling between weakly interacting Bose-Einstein condensates 
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Based on a tunneling Hamiltonian description, we calculate the Josephson, normal and interference 
currents between two Bose-Einstein condensates described by the Bogoliubov theory. The dominant 
Josephson term is of first order in the tunneling with a critical current density proportional to 
the ground state pressure. In contrast to superconductors, the normal current remains finite at 
zero temperature. We discuss the dynamics of the relative phase in a semiclassical approximation 
derived from an exact functional integral approach, which includes the interaction effects at fixed 
total particle number, ft is shown that the normal current leads to a damping of the Josephson 
oscillations and, at long times, eliminates the macroscopic quantum self trapping predicted by Smerzi 
et.al. Finally we give estimates for an experimental realization of Josephson tunneling in cold atomic 
gases, which indicate that coherent transfer of atoms might be realized with a 23 Na condensate. 
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I. INTRODUCTION 

The realization of Bose-Einstein condensation (BEC) 
in atomic vapors provides an example of superflu- 
idity, to which the approximation of a weakly interacting 
Bose gas applies quantitatively. Due to the low densi- 
ties n « 10 14 cm -3 and with typical scattering lengths 
a w 5nm, the gas parameter na 3 is of order 10 -5 , i.e. well 
in the range of applicability of the Bogoliubov theory. 
Although the interactions still entail strong quantitative 
changes compared to an ideal gas picture [Q , the systems 
are rather well described by the weak coupling Gross- 
Pitaevski equation (GPE) and its small fluctuations, the 
Bogoliubov equations [^). Many of the phenomena pre- 
dicted within this framework like sound modes with a 
linear spectrum have been established experimentally || . 
More recently genuine superfluid properties like quan- 
tized vortices [f7|8| or the existence of a critical velocity 
H have also been observed. 

The Josephson effect is one of the prime exam- 

ples of macroscopic quantum effects, displaying directly 
the broken symmetry associated with the relative phase 
of two weakly coupled condensates jL2fl . In supercon- 
ductors it is a rather standard phenomenon, in contrast 
to Bose condensed systems, where it has first been ob- 
served by Avenel and Varoquaux in superfluid 4 He [ jT3| . 
For dilute atomic gases, coherent oscillations have been 
seen in driven two component BECs |j~4| and in verti- 
cal arrays of traps of an optical lattice However 
these are Rabi- or Wannier-Stark-type oscillations and 
should not be confused with a genuine ac-Josephson ef- 
fect as discussed here. Theoretically, the Josephson effect 
in this context has been investigated by Smerzi et al. |l6[ , 
who focussed on the nonlinear dynamics of the relative 
phase and population difference. Zapata et al. |l7j calcu- 
lated the Josephson coupling energy due to condensate- 
condensate tunneling for weakly coupled harmonic wells 
and gave estimates for the normal currents due to the 



thermal cloud. These studies were essentially limited to 
a mean field like Gross-Pitaevski description, neglecting 
noncondensate contributions. More recently Villain and 
Lewenstein |0| showed that particles out of the conden- 
sate lead to a damping of the relative phase in a two 
component condensate with a Josephson like coupling in- 
duced by the external driving field juj . 

Our aim in this work is to provide a microscopic cal- 
culation of the Josephson effect between two weakly in- 
teracting BECs separated by a tunnel barrier (l9). Us- 
ing an idealized model, we give quantitative results for 
both the superfluid and normal currents per area which 
only depend on the microscopic parameters of the con- 
densate and the barrier transmission amplitude. More- 
over, a complete theory is given for the interaction effects 
arising in coupled condensates with a fixed total parti- 
cle number in the limit, where the number of exchanged 
particles is small. The associated dynamics is discussed 
within a semiclassical approximation. Using realistic pa- 
rameters in dilute atomic gases, the effects due to a fixed 
number of particles are appreciable, but do not destroy 
the qualitative structure of oscillating Josephson currents 
for condensates at different chemical potentials. Our es- 
timates indicate that the Josephson effect might be ob- 
servable with currently available techniques for 23 Na or 
lighter atoms. 

The outline of the paper is as follows: In section 2, we 
introduce the tunneling Hamiltonian model, including a 
discussion of the relevant states in condensates with a 
fixed total particle number and the general form of the 
Josephson coupling energy. This energy and the asso- 
ciated nondissipative current is calculated up to second 
order in the tunneling matrix elements in section 3. For 
a finite, given difference in the chemical potentials we 
determine in section 4 the complete particle current in 
terms of the normal and anomalous spectral functions of 
the individual condensates. In section 5, we develop an 
exact functional integral representation of the reduced 
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quantum mechanics of weakly coupled Bose condensates 
which incorporates the effect of interactions at fixed to- 
tal particle number. In section 6, explicit results are 
given for an idealized geometry, determining the physi- 
cally relevant currents per area from microscopic param- 
eters which characterize the separate condensates. The 
resulting semiclassical dynamics of the relative phase, in- 
cluding 'charging' effects and dissipation, is solved nu- 
merically in section 7 for parameters, which are realistic 
for current condensates of dilute atomic gases. Moreover, 
we test the tunneling Hamiltonian Ansatz by comparing 
it with exact numerical results for the one-dimensional 
GPE with a barrier. Conclusions and open questions are 
discussed in section 8. 



II. THE TUNNELING HAMILTONIAN 

Conceptually, the most obvious way to realize a 
Josephson geometry for BECs with current experimen- 
tal setups, is to split a single condensate in a long trap 
into two separate parts by a narrow light sheet produced 
by a blue detuned laser. The repulsive potential due to 
the ac-Stark shift is proportional to the laser intensity 
and thus the height and width of the barrier may be 
varied in a considerable range. In the limit of a strong 
barrier one has to zeroth order two completely separate 
condensates a and b with Hamiltonians -ff a /b- Provided 
that the coupling energy Ej due to the transfer of par- 
ticles across the barrier is small compared to the ground 
state energies of H a /bi which are of order fj,N, the cou- 
pled system may then be approximated by a tunneling 
or transfer Hamiltonian ppfl 

H = H a + H b + H T (1) 

in which the contribution Ht — A + describing the 
transfer of particles between a and b can be treated per- 
turbatively. Formally, introducing a set of orthonormal 
eigenstates {]/)} and {|r)} which decay exponentially in 
the barrier where a and b overlap, an instantaneous trans- 
fer is described by a term 

A = — ti r aja r (2) 
l,r 

with tunneling amplitudes ti r , which may be expressed 
in terms of current matrix elements between the left and 
right eigenstates p(J (see section 6 below) . The signs here 
have been choosen in such a way that for positive, real 
ground state wave functions in a and b, the associated 
tunneling amplitude is positive, giving rise to a lower- 
ing of the ground state energy in the coupled system, as 
expected in a standard double well situation. While a 
rigorous derivation of the transfer Hamiltonian is rather 
difficult [2(J , it is expected to be valid for high and narrow 
barriers, for which the mean field interaction of the con- 
densate is negligible within the barrier. This is in fact 



confirmed through numerical calculations in an exactly 
soluble one-dimensional geometry in section 7. Note that 
the states {\l}} and {|r)} introduced in Eq. (||) are not 
mutually orthogonal, and thus H a and H b do not com- 
mute beyond zeroth order in Ht 0]. Moreover, the 
creation and annihilation operators in (2) refer to the 
transfer of atoms, not that of quasiparticles, which are 
the proper excitations only within the individual conden- 
sates. 

As pointed out by Ferrell and Prange |2^] , the essence 
of the Josephson effect is that in the absence of any drop 
A/it in the chemical potential between both sides, no en- 
ergy is required to transfer a condensate atom across the 
barrier (note that a finite value of A[i may be present ei- 
ther eaternally e.g. by a drop in the gravitational poten- 
tial as in the experiment by Anderson and Kasevich [fL5| , 
or may arise internally at finite transfer currents through 
'charging' effects discussed in section 5 below). Thus the 
states \v) = \N a + v , N\, — v) in which an arbitrary in- 
teger number v of condensate atoms have tunneled from 
right to left, are degenerate. The ground state of the 
coupled system is therefore not an eigenstate \v) with a 
definite relative particle number, but rather a coherent 
superposition 

\ip)=Y,c v e^\v) (3) 

with a definite relative phase ip = ip a — tpi,. Here the 
c„ are real coefficients obeying ^2^ = 1 which are con- 
stant around v = 0, decaying to zero only at \v\ of order 
yJN a /b. In terms of the eigenstates \ip a , <fb) in which each 
of the condensates has a definite overall phase ip a .b, the 
state \tp) can be expressed in the form 

\V)« jf ^e-^+^ltpb + V^b), (4) 

using the standard relation between states with a fixed 
phase or a fixed particle number p2| . The state with 
a definite relative phase is thus a projection of states 
\cp a , iff,) with a broken gauge symmetry in the individual 
condensates to one with a definite total particle number 
N a + N b . 

The derealization in the relative particle number v 
described by the coherent state \<p), gives rise to a lower- 
ing of the energy of the coupled condensates in the same 
manner, in which a Bloch state gains kinetic energy by 
spreading over many sites in a periodic potential with lo- 
calized Wannier orbitals labeled by v. In a time reversal 
invariant situation, this energy (or free energy at T ^ 0) 
must be even and periodic in ip — > ip + 2n. Quite gen- 
erally it can thus be expanded in a Fourier series of the 
form H 

OC 

cos n<p. (5) 

n=0 

Since v and H<p are canonically conjugate variables, the 
nondissipative current 
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n— 1 

can be obtained from the knowledge of the phase depen- 
dent part of the coupling energy between the two con- 
densates. Within a tunneling Hamiltonian description, 
the coefficients E n = hl n /n are of order \Ht\ u ■ It is 
thus usually sufficient to keep only the lowest nonvan- 
ishing term, conventionally denoted by —Ej, which is 
negative in the standard situation, where the lowest en- 
ergy state is that with a vanishing relative phase ip = 0. 
As will be shown below, for Bose condensates the Joseph- 
son coupling arises already in first order in the tunneling 
amplitudes, 

< ip\H T \f >= -Ej cost/? (7) 

with a positive Josephson coupling energy Ej, which 
immediately determines the associated critical current 
I c = Ej/h. For superconductors, in turn, where the 
Josephson effect is associated with the coherent tunnel- 
ing of pairs, the first order term vanishes and Ej is pro- 
portional to the average of \ti r \ 2 at the Fermi energy. The 
latter being a direct measure also of the normal currents, 
one obtains the well known relation 

2Ef cs = A tanh (3A /2 • hG n /4e 2 (8) 

connecting the BCS Josephson coupling to the super- 
conducting gap A and the dimensionless normal state 
conductance G n 

Since tunneling is a small perturbation, the fact that 
there is a Josephson effect already in first order for a 
Bose condensed system, might seem to make higher order 
calculations unnecessary in this case. In fact, there are 
two reasons why this is not so: First of all, it turns out 
that realistic barriers which allow to observe a Josephson 
effect in coupled BECs have to be small enough, that 
higher order effects are nonnegligible. More importantly, 
though, it is only at second order that dissipative effects 
appear through normal currents which, as we will see, 
can have strong effects even if they are small. 

III. THE PHASE DEPENDENT COUPLING 
ENERGY 

In order to calculate the change in energy due to tun- 
neling, we employ perturbation theory up to second order 
in the transfer term Ht- To be specific, we assume that 
the condensate size is large compared with the coherence 
length £ = {Rima)" 1 / 2 over which the condensate wave 
function varies pj . Under this condition, which is well re- 
alized in present gaseous BECs, the eigenstates |Z) and \r) 
introduced above for a weakly interacting gas are given 
by 

• the condensate ground states <f> a /b(x/y) which fol- 
low from a solution of the Gross-Pitaevski equation 
in the given confining potential, and 



• a set of excited states orthogonal to (f> a /b( x /y)> 
which may be labeled by nonzero wavenumbers k or 
q for condensates a and b respectively. These states 
are obtained by solving the standard Bogoliubov- 
equations in a local density or - equivalently - a 
semiclassical approximation, resulting in wavefunc- 
tions and energies like that in a homogeneous sys- 
tem with a local condensate density n c (x) 

The transfer operator 

A = -t cc a\ a a^ h - t kc a\a c , b 
k 

- *cg4,a a 9 - t kq a\a q (9) 

q k,q 

is thus split into contributions describing condensate 
to condensate (c-c), condensate to noncondensate (c- 
nc) and noncondensate to noncondensate (nc-nc) tun- 
neling. To first order in Ht, only the c-c term con- 
tributes and gives rise to an energy gain of precisely 
the form (7), with a positive Josephson coupling energy 
Ef EC = 2t cc (N a N b ) 1 / 2 , which favors a fixed relative 
phase ip = in the ground state of the coupled system. 
Here, as in the rest of our work, we have assumed that \v\ 
remains much smaller than the average number N a / b of 
condensate atoms in each well. The Josephson coupling 
energy is thus independent of the number of transferred 
atoms to lowest order. 

For the calculation of the second order energy shift 

e 

involving all possible excited states |e), we follow a 
derivation given by Ferrell for the analogous case of 
Josephson tunneling between two BCS superconductors 
p6| . As is evident from (9), there are three contributions 
to the transfer Hamiltonian involving excited states: two 
c-nc terms and one nc-nc contribution. Denoting time 
reversed states by k = —k and using the time reversal 
invariance relation t ck = tj. c , the contribution which de- 
scribes tunneling between condensate b and nonconden- 
sate states in a can be written as 

fic b - nCa = _J2t kc (a{a cb + al b a k ) . (11) 

k 

Starting in a state \i>) with a definite number of atoms in 
each condensate but no excitations, the relevant excited 
states \k,v) associated with the first term in (11) are 
those, in which a condensate atom from b is transferred 
to a quasiparticle k in a. With the standard Bogoliubov 
transformation 

4 = M fc4 ~ v k & k ( 12 ) 

between the Boson or quasiparticle creation operators a) k 
or d k , the associated matrix element is 
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(k,v\H?- nc °\v) = -t kc u k ^W b , (13) 

again assuming \v\ <C N b . Now precisely the same excited 
state with energy E k (both E k and the particle and hole 
amplitudes u k and v k are positive and even in k, following 
standard notation |2?j) can be reached from the state 
\v + 2) by converting two condensate atoms in a to a 
pair (fc, k) and transferring one of the partners k to a 
condensate state in b, as indicated by the second term in 
(9). The amplitude for this process is 

(k, v\H c ^- nc «\ V + 2) = t kc v k y/W b . (14) 

It differs from (13) by the replacement u k — > — v k since 
the amplitude for creating a quasiparticle with momen- 
tum k by destroying an atom with momentum — k is — v k . 
In more physical terms, the relative minus sign between 
these amplitudes may be understood by noting, that in 
the Bogoliubov theory (and in fact quite generally for a 
Bose condensed system) there is a fixed relative phase 
7r between the condensate and pairs (fc, k) of nonconden- 
sate atoms p8|] . This minimizes the repulsive interaction, 
giving rise to a gapless excitation spectrum in contrast to 
the result of a Hartree-Fock approximation. The phase 
locking between condensate and noncondensate is also 
implicit in the form 

\<pa) cc ex P (e^y/K&U - e 2 ^ £ ^a£a^ (15) 

of the Bogoliubov ground state for a homogeneous sys- 
tem with a well defined overall phase tp a , as introduced 
in (4) above. Now, as emphasized before, the ground 
state of the coupled condensates is one with a definite 
relative phase <p, in which the amplitudes (13) and (14) 
add coherently 

(k, v\H c T »- nc « \tp) = -e iv n kc y/W b {u k - e 2 ^v k ) . (16) 

With E k as the relevant excitation energy, this gives rise 
to a phase dependent contribution 

^ b J2^-2u k v k -cos2ip (17) 

k k 

to the second order energy shift. Using 2u k v k = fi a /E k 
and exchanging the roles of condensates a and b, the 
total c-nc contribution to the phase dependent energy in 
second order is 




•cos2<^. (18) 



Obviously this term, which is much larger than the nc-nc 
contribution derived below, favors a relative phase ip = 
it/2 of the coupled condensates. In the weak tunneling 
regime considered here, however, it is always the leading 
first order term (7) which dominates, leading to a ground 



state with ip — 0. To obtain the nc-nc contribution, we 
use again time reversal invariance to write 

£nc-nc = _ £ ^ ( a t &g + a t % ) . (19) 
k,q 

The relevant excited states |fc, g, v) now have energy E k + 
E q and are characterized by N a +v and N b —i>—2 atoms in 
condensates a and b plus one quasiparticle k respectively 
q on each side. Corresponding to the two contributions 
in (19), these states can be reached starting from either 
\v) or \v + 2). Adding the respective amplitudes in the 
actual coherent ground state with a well defined relative 
phase, one obtains 

(*, q, v\H^- nc \ V ) = e^t kq (u k v q + e 2 ^v k u q ) . (20) 

The phase dependent part of the associated second order 
energy shift is thus 

A4t«c(^) = -2X; M 2 ^^.cos2^, (21) 

£j k + Ej q 
k,q H 

in perfect analogy to the result obtained for a super- 
conducting Josephson contact |24|,|26|], where this is the 
only contribution. Using the relation between the aver- 
age of \t kq \ 2 at the Fermi energy and the normal state 
conductance, (21) directly leads to the result (8) at zero 
temperature. For weakly interacting Bose gases, in turn, 
with the explicit results (53) and (54) for the matrix ele- 
ments given below, the contribution (21) is easily shown 
to be smaller than the c-nc contribution (18) by a factor 
Vna 3 |2{|, and thus is negligible for gaseous BECs. 

For a given static value ip of the phase difference be- 
tween both condensates, the nondissipative current ob- 
tained from (6) up to second order in Ht is thus of the 
form 

I(ip) = -I c sin (p - Ji(0) sin 2(p, (22) 

with a positive critical current I c = Ej/H. The magni- 
tude Ji(0) < of the second order contribution is — 2/ ft 
times the factor in parentheses in (18). The small second 
order term changes the value of ip at which the current 
is maximal from 7r/2 to ir/2 + 2| Ji(0)|// c , however the 
critical current itself is unchanged to this order. In an 
open system, with a reservoir of particles, where it is 
possible to impose a constant external current, the phase 
difference adjusts itself to a constant value, determined 
by (22). This is the well known dc-Josephson effect, i.e. 
a finite current at vanishing chemical potential difference 
A/i = 0. In the case of two coupled BECs with a fixed 
total number of particles, any current flow is connected 
with a finite value of A/i even in the absence of an ex- 
ternal potential drop, because [i depends on the particle 
number. Thus, by the Josephson relation (23), one has 
inevitably a phase difference which evolves in time. In 
addition, finite dissipative currents appear, requiring a 
fully dynamical treatment, as will be given in the follow- 
ing sections. 
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IV. JOSEPHSON AND NORMAL CURRENTS AT 
GIVEN Afj, 



The first term on the right hand side of Eq. ( fjlf ) is 
then 



In the following, we want to determine the complete 
current in a situation with a finite difference in the chemi- 
cal potentials. We start by considering the idealized case 
in which A/z is considered as a fixed, externally given 
value. This approximation, which applies to open sys- 
tems like standard Josephson contacts between super- 
conductors, neglects the change in the chemical potential 
associated with the transfer of particles in a system with 
a fixed total number of particles, an effect which is taken 
up in section 5. A finite value of A/z gives rise to an ad- 
ditional term A/j, v to the total Hamiltonian (1). In close 
analogy to the calculation of currents in normal and su- 
perconducting tunnel junctions, this term may formally 
be eliminated by a time dependent gauge transfomation. 
One thus obtains a tunneling Hamiltonian with A/i = 0, 
in which the transfer matrix elements are modulated in 
time with a factor e llp ^ > such that 



*■ d , , 



-A/i. 



(23) 



Since the current operator 

I=~N a = i{k-lS)/h (24) 

is already linear in the tunneling matrix elements, the 
problem can, up to second order in Ht, be formally 
treated like in time dependent linear response pofl . It is 
then straightforward to show that the expectation value 
I of the current is given by 

j = I{ e -W')A(t)-e i * , WA t (i)) 
h 

+ / dt'{e- i( *'W- (p(t ' )) ([A(t),A t (« , )]> (25) 



-i((p(t)+tp(t')) 



<[A(t),A(f)]» 



where the time dependence of the operators A and has 
to be taken with respect to the unperturbed Hamiltonian 
H a + Ht. The expectation values reduce to ones in the 
associated ground state, provided we restrict ourselves to 
the limit of zero temperature, as is done throughout in 
the following. 

As was pointed out above, the eigenstates of H a + 
for gaseous BECs can be choosen as the Gross-Pitaevski 
wave functions in each well, and a set of excited states 
obtained from solving the Bogoliubov equations in local 
density approximation. As long as the number of trans- 
ferred particles remains small compared to N a / b , these 
states are unaffected to lowest order. Thus both the con- 
densate wave function as well as the quasiparticle en- 
ergies E k and amplitudes u k and v k are unchanged by 
the current flow. The comparison with the numerical 
results in section 7 will show, that at least on the Gross- 
Pitaevski level this approximation is well justified in the 
weak tunneling limit discussed here. 
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-<*>(*) A(t) -e^W A+ (t)) 



2t r 



N£N§ sin ip{t) = -I c sin tp(t) , (26) 



in perfect agreement with the result derived in the pre- 
vious section for a static situation. The second term in 
Eq. ( p5|) involves noncondensate tunneling. At fixed A/z, 
it is readily evaluated within Bogoliubov theory, giving 
three contributions: 



h 2 



di'{ e -^(*)-?(0)([A(i),At(/)]) 

+e -*M*)+¥>(0)([A(t),A(t')])} 



= J„(A/i) - Ji(A/x) sin2<p(t) + J 2 (A/z) cos2^(i). (27) 

The first one is a phase independent normal current, ac- 
counting for the expected dissipative current flow towards 
the lower chemical potential. The sin 2(/?-term describes a 
second order Josephson current, associated with the cor- 
responding phase dependent coupling energy derived in 
the previous section. Finally there is an interference cur- 
rent proportional to cos2y> which is out of phase by n/2 
from the corresponding Josephson contribution. Similar 
to the related cos 99-term in superconductors |ll| , this is 
a dissipative current which vanishes at A/i = 0, as the 
normal current does. Within the local density approx- 
imation for solving the Bogoliubov-equations, the cur- 
rent amplitudes J n , J\ and J 2 can be expressed in terms 
of the tunneling matrix elements and the normal and 
anomalous spectral functions A(k,uj) and B(k,u>) of the 
weakly interacting homogeneous Bose gas. In standard 
notation |27| they read 

A(k, to) = 2n{ul8(w - E k /h) - v 2 k 6(uj + E k /h)} (28) 



B(k,u) = 2iru k v k {5{tj + E k /h) - 5{lu - E k /h)} . (29) 

With ns(x) = (exp/3x — I) -1 as the Bose distribution 
function, the current amplitudes in Eq. ( p7| ) are then 
given by 



J„(A M )=A a c £ M^A/Vft) 



k 

A(k, uj)A(q,uj + A/j,/h) 



\t kc \ 2 

h 2 



A(k,-Afi/h) 



(30) 



and 
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J 2 (A/x) = N%J2 l -^B(q,-An/K) 



k 

I^nfJZ nB{hw + A/i)] 



B(k,u))B(q,u + Afi/H). 



(31) 



The amplitude Ji(A/i) of the second order Josephson 
current can be obtained from that of the associated dis- 
sipative contribution J 2 (A/i) via a Kramers-Kronig rela- 
tion 



Ji(Afi) = -V 



duj J2{huj) 

7T LU — A/i/ft 



(32) 



It is straightforward to check, that at zero temperature 
and vanishing A/i, the combination hJ\(0)/2 is identical 
with the sum of the coefficients of cos 2<p in (18) and (21) 
of the second order phase dependent coupling energy, as 
it should. As noted in the previous section, a particular 
feature of tunneling between Bose condensates is the ap- 
pearance of c-nc contributions. Indeed, as the first two 
terms in Eq. ( |30| ) show, they also contribute to the phase 
independent normal current, where they are in fact dom- 
inant compared to the conventional nc-nc contribution 
described by the last term. Another point that should 
be mentioned, is the absence of a second order contri- 
bution due to c-c tunneling [p9[ . This follows because 
the commutators involving only condensate operators in 
Eq. (|2^) vanish. A different proof of this result may be 
obtained, by considering all possible Feynman diagrams 
associated with second order tunneling processes. Dia- 
grams proportional to \t cc \ 2 are reducible and thus give 
no separate contribution. 

The similarity of the second order currents in BECs 
and superconductors is a result of the close formal anal- 
ogy of Bogoliubov and BCS theory. The additional con- 
tributions due to c-nc tunneling in BECs lead to inter- 
esting new features, however. As will be discussed in sec- 
tion 6, for the Bose case the dissipative currents J„(A/i) 
and J 2 (A/x) remain finite at zero temperature in contrast 
to superconductors, where they vanish exponentially due 
to the presence of an energy gap. These currents thus 
provide a mechanism for the damping of Josephson os- 
cillations even at T = 0. Similar to the situation in 
superconductors ]3C| ] , it may be shown that the total dis- 
sipation associated with the two irreversible currents is 
always positive p9| ]. Physically, this property is equiv- 
alent to the statement, that the sum J n + J 2 cos 2ip is 
strictly positive for A/x > 0, i.e. the total dissipative 
current always flows towards lower chemical potential. It 
therefore leads to a reduction of A/i and hence to a damp- 
ing of the phase dynamics, similar to the contribution of 
excited states in driven two component condensates dis- 
cussed in 18 1 . 



V. PATH INTEGRAL DESCRIPTION 

So far, we assumed hip to be given externally, thus ne- 
glecting the variation of the chemical potential with the 
total number of particles. As we will see, however, this ef- 
fect is quantitatively important even in the regime where 
N a — Nf, <C N a /b- To lowest order in v/N a /b, the depen- 
dence of the chemical potential on the actual number of 
particles can be accounted for by an additional 'charging' 
energy U(N a - N k ) 2 /8, with U = dfi a /dN a + dfx b /dN b 
being of order fi/N. In the limit U > HI C , this term sup- 
presses phase coherence, since by canonical quantization 
v — > id I dip , charging effects are equivalent to quantum 
fluctuations of ip. For any U ^ 0, the chemical poten- 
tial difference in hip = —A/i contains an internal con- 
tribution Uv proportional to the number of transferred 
bosons. The resulting dynamics of ip is therefore intrinsi- 
cally nonlinear. Within a simple description, these effects 
have already been introduced by Smerzi et.al. Q, who 
also took into account the (/-dependence of the Joseph- 
son coupling Ej cx \J (N a + v) (Nb — v) . While this effect 
is indeed relevant for strong asymmetries v = 0(N a /b), 
it is negligible compared to the charging term of order 
\iv 2 I 'N provided we are in the regime Ej <C iiN where 
Josephson tunneling is a small perturbation. In order to 
include charging effects in a microscopic description, we 
use the formalism developed by Ambegaokar, Eckern and 
Schon fjlf for superconductor tunnel junctions. 

In a coherent state path integral formulation for inter- 
acting bosons, the grand- canonical partition function of 
two coupled condensates is given by 



DtpaDtpb exp - 



0h dr 

„ h 



d 3 xipl(hd T - fJ, a )lp a 



(33) 



Here [i a jb fix the average number N a /b of atoms in each 
condensate, while H(ip a ,il>b) is the standard representa- 
tion of the basic Hamiltonian (1), in which the bosonic 
field operators are replaced by complex c-number fields 
4>a/b(x/y, T ), depending both on the spatial coordinates 
plus an imaginary time like variable r [ B2| . With the 
charging term explicitly included in H(ijj a ,ipb), Eq. ( [33] ) 
is also appropriate for describing a system with N a + Nb 
fixed. As in the analogous fermionic problem of coupled 
superconductors, it is convenient to introduce an aux- 
iliary path integral J DV to remove the charging term 
quartic in the fields via a Hubbard-Stratonovich trans- 
formation. The functional integral is then evaluated in a 
saddle point approximation, the saddle point being just 
the condensate wavefunction for two separate conden- 
sates. Performing a gauge transformation, we may re- 
place the integration variable V(t) by the phase ip(r), 
which is the dynamical variable of interest. Further- 
more, we introduce Nambu spinors = (ip*, ip a) i/j^, tpb) 
for the Gaussian fluctuations around the saddle point 



G 



<J>^ = (tfia,(f) a ,(j)l,(j)b). The full partition function is then 
given by 



Z = e 



exp - 



Dip J d^d¥ 
7»V 



(It 



2U 



Ej cos ip 



+- J d s xd s y¥{-Cr l + i)# + #t£$ + $t£j> 

(34) 

where So is a bulk contribution to the action. In coordi- 
nate space, the 4x4 tunneling matrix t is defined by 



t = 



t X y 



t = 



t xy 
t 



xy 



while the differential operator G satisfies 
-l ( G a \ _ { 5(x- x')l 



G~ 



G 



6(y - y')l 



(35) 

8{t-t'). 
(36) 



In the last equation, G Q/ / h is shorthand for the noncon- 
densate part of the 2x2 matrix Green's function for a 
weakly interacting Bose g defined in [E7|. 

Completing the square in the exponent of (|34|), the 
Gaussian integral over \I/ is readily performed and yields 



Z a Z b / D(p(T)exp(-S[(p]/h). 



(37) 



Here, Z a j b are the partition functions of the separate con- 
densates, while tunneling effects are taken into account 
by an effective action for the relative phase 



S[ V ]: 



-h I dr dr' [a(r — t') cos(tp — tp') 



o 



+/3(t - t') cos(v? + ip')] 



(38) 



where tp — ip(r) and <p' = V 5 ( r ')- The integral kernels 
a(r) and /3(r) may be expressed in terms of the currents 
J n (A(i) and J2(A/i), introduced in the previous section 
via their Fourier coefficients 



(3(iu n ) 



duj J n (—hu>) 
2n iuj n — lo 
duj J2(huj) 
2tt iu) n — w ' 



in an expansion of the /3?i-periodic functions 



(39) 
(40) 

(41) 



in bosonic Matsubara frequencies uj n = 27rn//3h, n inte- 
ger. 

The effective action Eq. (Q) is formally identical with 
the one obtained for a superconducting Joscphson con- 
tact |3^|, except for the additional local contribution 
— Ej cos <p{r) describing the first order Josephson effect 
due to c-c tunneling. The terms nonlocal in time give rise 
to both the second order Josephson currents related to 
P(t—t') and the normal current associated with a(r— r'). 
For the idealized model discussed in section 6, where 
Ja(A/i) and J„(A/x) are linear for small values of the 
chemical potential difference, both a(r) and /3(t) asymp- 
totically decay like 1/r 2 , i.e., the Josephson contact ex- 
hibits ohmic dissipation. 

An exact evaluation of the path integral in (37) is obvi- 
ously impossible, since the effective action is nonlocal and 
- in particular - nongaussian. Nevertheless, as shown by 
Ambegaokar, Eckern and Schon pl| , it is a useful start- 
ing point for the derivation of a semiclassical approxima- 
tion for the phase dynamics. Indeed, from the imaginary 
time description, it is possible to derive a semiclassical 
equation of motion in real time, which has the form of a 
classical Langevin equation with state dependent quan- 
tum noise. In the limit where J2(A/x) and the normal 
current J„(A/i) = G„A/i are purely ohmic, it reads 

-^j- + G(ip)hip + I c sin ip + Ji sin 2tp = rji cos ip + 772 sin <p, 

(42) 

where 771.2 are independent Gaussian random forces. 
Their autocorrelation function may be expressed in terms 
of the real time functions a(t) and f3(t) via 

<r?i/2(*)r?i/2(*')) = 2[« 7 (* - t') T P\t - t')] , (43) 

where a 1 and f} 1 are the imaginary parts of the analytic 
continuation of a(r) and /3(r) similar to plf . The damp- 
ing term is proportional to a phase dependent effective 
conductance G(<p) = G n (l + cos 2<p), with a normal con- 
ductance G n , which will be explicitly evaluated in section 
6. It gives rise to a relaxational dynamics with a typical 
time scale l/G n U. Even without damping, however, the 
tp-term leads to currents which are not perfectly sinu- 
soidal for a given external chemical potential difference. 
Neglecting the fluctuating forces, the resulting average 
dynamics will be discussed in section 7 for realistic val- 
ues of the parameters in weakly coupled BECs. 



VI. AN EXPLICIT MODEL 

If the trap potential close to the barrier varies slowly on 
the scale of the coherence length £, we may apply the ap- 
proximations discussed in section 3, using eigenstates of 
a locally homogeneous system with constant external po- 
tential throughout a and b, respectively. A simple model 
geometry for such a Josephson junction is shown in Fig. |]: 
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FIG. 1. The model geometry considered for the evaluation 
of the current densities. 

Two BECs confined in a cubic volume L 3 are separated 
by a square potential barrier of height Vb and width d. 
For simplicity, we impose periodic boundary conditions 
in the directions parallel to the potential barrier. This 
simple model accounts for the fact, that all currents must 
scale linearly with area. It allows us to calculate the rel- 
evant currents per area in terms of microscopic parame- 
ters of the bulk condensates. For a concrete experimental 
setup, one may then determine the actual critical current 
I c = j c A from the associated current densities and the 
effective contact area A. 

For this system, the currents I Cl J„, J\ and J2 may 
be evaluated explicitly from the detailed form of the tun- 
neling matrix elements. Since we have assumed trans- 
lation invariance parallel to the barrier, the associated 
momenta fen , q\\ are conserved, and the problem is effec- 
tively one-dimensional (Id). The tunneling amplitudes 
are then obtained from the current matrix elements [EOl 



2m 



XV 



dXr 
dz 



dxi 
dz 



(44) 



of the Id wave functions \^ r (z), taken at the center of 
the barrier at z = 0. In the limit of weak tunneling, 
the barrier height Vb is much larger than the average 
chemical potential p, — (fi a + ^b)/2. The wave functions 
in (44) are therefore effectively single particle eigenstates 
which decay exponentially like expift^z, with inverse 
characteristic length = y // 2m(V B — fJ-)/fi 2 . For a small 
difference A/x = fi a — j-ib <C \x in the chemical potentials, 
(44) then reduces to 



tlr 



h Kr, 



(45) 



The calculation of the different tunneling amplitudes t cc , 
tk c and tk q , thus requires 

• a solution of the Id GPE for a finite barrier, which 
approaches 4> c — \/N c ■ \c —> \fn~c far from the 
barrier, i.e. a uniform condensate with density n c 
(note that the boundaries at z = ±(L + d/2) are 
irrelevant since L>^), and 



• corresponding solutions of the Id Bogoliubov equa- 
tions with finite momenta k , q, which smoothly 
connect to the exponentially decaying single par- 
ticle states below the barrier. 

For concreteness we consider the eigenstates \i m con- 
densate a, to the left of the barrier. It is then straight- 
forward to show that 



1 e 



-K Ma (z+d/2) 



-K t , a (z+d/2) 



(46) 



provides an appropriate solution of the GPE. It 
is obtained by connecting the well known solution 
-y/n^tanh ^| of the full GPE in the regime z < -d/2 to 
an exponentially decaying solution of the linearized equa- 
tion below the barrier z > —d/2. The condition that x 
and its first derivative are continuous at z — —d/2 fixes 
the prefactor in the solution of the linear equation. The 
linearization is justified in the limit k£ 3> 1, where the 
mean field interaction is negligible compared to the re- 
pulsive potential separating the two condensates. Using 
the solution (46) and the condition A/i <§; /i, the c-c 
tunneling amplitude is given by 



t cc = f(V B /ft)-^e^ 



KfiL 



with a correction factor 



f(x) = (l-(x-Vx 2 -l)) 2 , 



(47) 



(48) 



which is smaller than one and approaches unity in the 
high barrier limit Vb S> ft- Since t cc is proportional 
to l/L, the resulting Josephson coupling energy Ej — 
2t cc \/N a Wb scales like the contact area A = L 2 , as it 
should. The associated critical current density is 



j c = 2f{V B /n) z e 

tiKn 



(49) 



As expected, it is linear in the barrier transmission ampli- 
tude exp— Kpd. More interesting is the prefactor which, 
by the standard relation 2p = [in c for the ground state 
pressure of a weakly interacting Bose gas, is just pro- 
portional to the geometric average of the ground state 
pressures p a and pb of the two condensates. It is in- 
structive to compare this with the corresponding result 
for coupled ideal Bose gases, where the condensate wave 
functions are the well known wavefunctions of the one 
particle ground state in a potential well of finite height. 
The associated c-c tunneling amplitude is then readily 
evaluated, giving 



e) = 2^ 



1 



TO KpL 3 



(50) 



8 



^From this, we immediately find a critical current 



J<°> = 4tt 



(51) 



which is independent of the contact area. Evidently, the 
behaviour implied by Eq. (^lj) is an artefact of the ideal 
Bose gas, which exhibits a vanishing quantum pressure 
below the condensation temperature T c in the thermo- 
dynamic limit. The critical current density thus remains 
finite for large contact areas only if the repulsive interac- 
tion between the particles is included. The linear scaling 
of the critical current with area is also found in the model 
by Zapata et al. |L7|]. In fact, our result for the critical 
current agrees with theirs for a plane barrier to leading 
order in e~ K > id , i.e. with / = 1. As a final point of our 
discussion of c-c tunneling, we mention that due to the N- 
dependence of the chemical potential in (47), the ampli- 
tude t cc - which is identical with the coupling parameter 
K introduced in |l6f| - depends on the number of particles 
in the condensate. In the limit \v\ -C N a jb discussed here, 
this effect is negligible compared to the charging effects, 
as mentioned above. In situations, however, where \v\ 
becomes itself of order N a / b , this dependence is impor- 
tant and the c-c tunneling amplitude cannot be treated 
as a constant, as assumed in ref. |Tg[ . 

To calculate the remaining amplitudes t^ c and tkq 
which enter into the higher order and dissipative cur- 
rents, we need to match a solution of the Id Bogoliubov 
equation for a homogeneous condensate in z < —d/2 to 
an exponentially decaying wave function below the bar- 
rier. Since the condensate density n c (z) there vanishes 
very quickly like exp— 2kz, the mixing between the am- 
plitudes Uk{z) and Vk(z) may be neglected (for simplicity, 
k and q here are just the z-components of the original 3d 
wave vectors k and q). Connecting the trivial solution 

1/2 

(2/L) ' sin (kz + S) for Uk(z) in a, to one below the bar- 
rier such that the function and its first derivative are 
continuous at z = —d/2, we find 



9\ 1/2 k 



(52) 



Here, it has been assumed that k <C K, i.e. the kinetic 
energy for motion perpendicular to the barrier is small 
compared to Vb- The tunneling amplitudes involving 
noncondensate states now follow easily from (45) and are 
given by 

h 2 k 

the = —z 7 e~^ d • 6t * (53) 

mtaKpL ^1 + (jfe£)2 *!I.O V ' 



and 



2h 2 



kq 



-_ e - K * d ■ St , (54) 



again assuming A/i <C fi. These matrix elements vanish 
linearly with the incoming momentum at low energies, as 



expected for a tunneling amplitude. In order to avoid the 
related unlimited increase of the amplitudes for large mo- 
menta, we have introduced a cutoff 1 / y/l + (k^) 2 which 
leads to a saturation of the matrix elements at an en- 
ergy scale still much smaller than the barrier height Vb- 
In fact for energies in this range, the transfer Hamilto- 
nian fails, because the coupling between the two sides 
can no longer be treated perturbatively. Fortunately, for 
small chemical potential differences A/k <C p, only the 
low energy phonon like excitations with sound velocity 
c = yj fi/m contribute to J n and J 2 - Indeed, for small 
frequencies, the spectral functions are antisymmetric in 
u) and are equal up to a sign 

A{k,u)~-B(k,u)='^5{uj-c a k) (w>0). (55) 

The sum over momenta is thus restricted to a regime 
where the cutoff is irrelevant. It plays a role only for the 
less important second order Josephson current J\ , which 
is small compared to the leading first order term, with a 
typical magnitude 



J 1 (A A i = 0)«/ c -e"^ rf /^- 



(56) 



Using (55), the c-nc contributions to J„ and J2 are read- 
ily evaluated in the limit L — > 00, where the one particle 
energy spectrum becomes continuous. To lowest order in 
the chemical potential difference, both currents are linear 
in A/i, behaving like 



MAfj. -»■ 0) = J 2 (A[i -f 0) = GnAfj,. 



(57) 



The associated normal conductance per area for a sym- 
metric situation is finite at zero temperature and given 
by 



G n /A 



2V2n c 



(58) 



Concerning the nc-nc terms, it is straightforward to see 
that their contribution to G n vanishes like T 4 and thus 
is negligible. The same applies to the nc-nc contribution 
to Ji(A/i = 0) which, as noted in section 3, turns out to 
be smaller than the c-nc term (56) by a factor V na 3 [ p9| . 

A remarkable feature of tunneling between Bose con- 
densates is the presence of a normal conductance pro- 
portional to the condensate density. Contrary to the sit- 
uation encountered in superconductors, the dissipation 
in a Bose Josephson junction thus remains finite even at 
T = 0. It is obvious that this result relies on the con- 
tinous excitation spectrum of the homogeneous system, 
and thus its relevance to BECs in traps might be ques- 
tioned. I practice, though, for Afi > tuo and for the 
experimentally relevant temperatures where ksT 3> fuj, 
the discreteness of the spectrum in the harmonic traps 
is irrelevant, and the continuum approximation thus well 
justified. We shall see, in the next section, that dissipa- 
tive currents play an important role in the dynamics of 
coupled BECs, even if they are small. 







VII. DYNAMICS FOR REALISTIC 
PARAMETERS 

On the basis of our microscopic model developed in 
the previous sections, we now aim to provide realistic es- 
timates for the necessary requirements and the expected 
dynamics in a Bose Josephson junction, which may be re- 
alized experimentally with present techniques. Assuming 
a contact area A = 20 /tm 2 and a condensate extending 
about 10 /im in the direction perpendicular to the bar- 
rier, we have roughly N a = Nb = 10 5 atoms on each side 
on average. With the known mass and scattering length 
the associated mean density b = 5x 10 14 cm -3 then de- 
termines the average chemical potential /2 = h x 8.2 kHz 
for 23 Na, which is the more favorable case for observ- 
ing tunneling due its lighter mass compared with 87 Rb. 
Since the tunneling amplitude decreases exponentially 
with barrier thickness d, it is desirable to have widths 
d as small as possible. With barriers realized by a blue 
detuned laser, d = 1 /jm is difficult to reach, but still a 
realistic value. Using (49), the resulting critical current 
can now be determined for any given barrier height. With 
a large barrier Vg = 5/2, however, the resulting critical 
current is only I c = 584 sec -1 , which is far too small to 
be observable. Indeed, the corresponding tiny Josephson 
coupling energy is not sufficient to establish phase co- 
herence across the barrier against even the minute ther- 
mal fluctuations at the typical temperatures of gaseous 
BECs. In fact, these fluctuations give rise to a smearing 
(A(p 2 ) = ksT/Ej of the relative phase. A Josephson ef- 
fect, which requires small phase fluctuations (Aip 2 ) -C 1, 
is thus only possible if I c ^> fcgT/fi, ~ 10 4 sec -1 at 
T = 100 nK. For conventional superconductors, where 
T is of the order of 1 K, this condition expresses the well 
known fact that minimum critical currents are around 

1/iAli- 

In order to observe Josephson tunneling in gaseous 
BECs, one therefore needs barriers which are not much 
larger than the average chemical potential. For a quan- 
titative estimate, it is important to consider the neces- 
sary magnitude of the oscillating currents, which are de- 
tectable as a clear signal of the Josephson effect against 
experimental noise. For a given difference A/i in the 
chemical potentials, the amplitude of the oscillation in 
the number of particles is Ej/Afi. With current ex- 
perimental resolution, this amplitude should be at least 
around five percent of the total particle number to be 
detectable. Using typical values A/i = 0.02/2, this re- 
quires Josephson coupling energies Ej > 10 -3 (N a + Nb)p. 
(note that this is still consistent with the requirement 
Ej <C (N a + Nf,)p, of a tunneling Hamiltonian descrip- 
tion). With the numbers above, this translates into criti- 
cal currents I c > 10 7 sec -1 , putting a much stronger limit 
than that of negligible thermal fluctuations. To realize 
a Josephson contact in the required range, it is neces- 
sary to choose a rather small barrier height Vb — 1.25/2. 
The critical current in a condensate with the above pa- 



rameters is then around 4 • 10 6 sec -1 . For 23 Na, the ob- 
servation of a Josephson effect is therefore close to the 
limit of detectability, and may be possible in an opti- 
mized setup. By contrast, for 87 Rb, the achievable criti- 
cal currents even for barriers as low as 1.25/2 are an order 
of magnitude smaller than those in 23 Na. Thus it appears 
unlikely, that coherent transfer of atoms across a tunnel 
barrer may be realized in this case with presently avail- 
able condensates. 

For a quantitative analysis, it is necessary to con- 
sider the time evolution of the relative particle number 
v = (N a — Nb)/2 including 'charging' effects and dissi- 
pation. As noted in section 5, the total difference A/t 
in the chemical potential is a sum of a possible external 
contribution A/to (induced e.g. by a difference in the 
gravitational potentials as realized in p5[ ) plus an in- 
ternal 'charging' term U v associated with a difference in 
the densities due to particle transfer. The generally valid 
Josephson relation (23) thus reads 

hip = -A/to ~ Vv . (59) 

For a constant A/to, a further time derivative then leads 
to the semiclassical equation of motion 

^j(p{t) +G n htp(t)(l +cos2^(i)) 

+I c sinip(t) + J 1 (-h(p{t))sin2ip(t) = (60) 

provided — z> is replaced by the total current obtained 
for a given instantaneous A/i = — tup. Obviously this is 
just Eq.(42) without the fluctuating forces, determining 
the dynamcis of the averaged relative phase and the as- 
sociated mean current. Eq.(60) is the analog of the well 
known resistively shunted junction model for supercon- 
ducting Josephson junctions p3|], and is valid for small 
phase fluctuations (Aip 2 ) <C 1. Apart from the con- 
dition Ej ^> ksT mentioned above, this also requires 
that the Josephson coupling is much larger than the 
charging energy U. Now, in a symmetric situation with 
N a = Nb — N and for the parameters choosen above, 
U = 2/2/jV is only of order h x 0.16 Hz compared with 
Ej = h x 0.6 MHz. Nevertheless, charging effects are far 
from negligible, since for small values of A/to they are the 
dominant contribution in (59) driving the phase dynam- 
ics. To give a quantitative example, we have numerically 
integrated the averaged, semiclassical time evolution of tp 
and v for the above parameters, neglecting second order 
currents. As is shown in Fig. 2, the normalized particle 
number difference z — vjN for a constant offset poten- 
tial A/to = 0.03/2 and initial condition t/(0) = exhibits 
appreciable deviations from a pure sinusoidal oscillation 
predicted in the absence of 'charging' effects. Obviously 
both the amplitude and the frequency of the Josephson 
oscillations are modified, although only by a factor of or- 
der of unity. The typical time scale for the oscillations 
is msec, which is well within the reach of experimental 
observation. 
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FIG. 2. Nonlinearity of the Josephson oscillations due to 
charging effects. For the parameters discussed in the text, the 
normalized population difference is shown for initial relative 
phases ip = 0, 1, 2 (from top) both with (solid line) or without 
charging effects (dashed line). The time is measured in msec. 

Regarding the higher order and dissipative currents, 
which have been neglected in Fig. 2, it turns out that for 
the low barriers choosen here, the amplitude of the sec- 
ond order Josephson current is not much smaller than the 
leading term, with a typical value Ji(A/i = 0) ~ 0.1I C . 
Now although this contribution oscillates at twice the 
frequency of the first order term, an experimental detec- 
tion is probably out of reach at present, given the narrow 
margin for seeing even the dominant effect, as discussed 
above. The dissipative normal and interference currents, 
however, significantly change the dynamics of the cou- 
pled condensates and thus should be indirectly observ- 
able even if the individual contributions are tiny. Indeed 
from Eq. (^50|) it is evident, that these currents lead to a 
damping of the phase dynamics with a typical time scale 
of order r = 1/G n ll. For the parameters choosen above, 
the normal conductance is G n = 250//i, giving a relax- 
ation time r = 0.025s. The amplitude of the Josephson 
oscillations will therefore decay to 1/e of their inititial 
value after typically 25 periods. As pointed out by Ru- 
ostekoski and Walls [M , dissipative currents also elim- 
inate the macroscopic quantum self trapping predicted 
by Smerzi et.al. |l6| , |35|| for systems with a large initial 
population imbalance and small Josephson coupling. For 
a quantitative estimate, we consider the influence of an 
ohmic normal current J n = G„A/i in a simplified equa- 
tion of motion for the normalized particle number differ- 
ence z. With t = 2t cc t/h as a dimensionless time and for 
a symmetric situation with A/iq = they read 



dz i 

= - z 2 sintp-2hG n fi/Ej ■ z (61) 

at 



and 



— = -Az- 
dt Z 



: COS if . 



(62) 



Here A = UN/2t cc = p,/t cc is essentially the ratio be- 
tween the condensate ground state energy E rs £LN and 
the Josephson coupling energy Ej = 2t cc N. The cosy- 
contribution in (62) arises from the z-dependence of the 
critical current I c (z) = I c • Vl — z 2 , relevant for large 
asymmetries z = O(l), while the second order Josephson 
currents discussed above are neglected for simplicity |36| . 
The phenomenological equations (61,62) lead to a self 
trapping of the condensate, provided A is larger than 
a critical value which depends on the initial asymmetry 
z(0). In the absence of dissipation, the asymmetry will 
be maintained in time, leading to a behaviour as shown 
in Fig. 3 (dotted line), with z(t) oscillating near its ini- 
tial value z(0) = 0.6 (as in ref. [IJ] we choose A = 11 and 

v(o) = o). 
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FIG. 3. Destruction of macroscopic quantum self trapping 
through dissipative currents. For comparison, the dynam- 
ics obtained without dissipation is also shown (dashed line). 
Time is measured in units of h/2t cc and thus with a typical 
value t cc — h x 20 Hz, t = 25 corresponds to about 0.1 sec 

It is obvious that normal currents, which lead to an 
equilibration of the chemical potentials A/i = fi a — fib — > 
for long times, will eventually destroy a self trapped 
state. Observation of macroscopic quantum self trapping 
is therefore possible only if the time scale l/G n U for equi- 
libration is much larger than H/2t cc , which is the typical 
scale for the dynamics in the trapped state. This requires 
that the conductance obeys hG n -C Ej/2ji, i.e. the co- 
efficient in front of the dissipative term in (61) should 
be small compared to one. Now from our microscopic 
results (49) and (58), we find that 



2hG n fi/Ej = 



2y/2 



(63) 



For high barriers, with a corresponding tiny Josephson 
coupling energy, this ratio can certainly be made small 
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enough, such that dissipation has essentially no influence 
on the dynamics. As pointed out above, however, realis- 
tic systems allowing to observe Josephson tunneling re- 
quire small barriers. For those, the inequality is strongly 
violated and in fact, for the parameters choosen above, 
the ratio in (63) is close to one. For a quantitative exam- 
ple, which allows a comparison with the results obtained 
in ref . fl6[| , we assume a condensate with much lower den- 
sity n m 1.3 • 10 13 cm -3 and a comparatively large value 
tec = h x 20 Hz, giving A = 11 [B7L With reasonable 
parameters Vb ~ 2/2, n^d — 4.2 and Kn£ — 3, the ratio 
in (63) has a rather small value 0.0275. The resulting 
time evolution of z(t) including dissipation, is shown as 
the solid line in Fig. 3. Evidently self trapping is de- 
stroyed rather quickly even in this case and in fact, the 
situation does not change qualitatively for other parame- 
ters which seem accessible. From our microscopic results, 
it thus appears unlikely that macroscopic quantum self 
trapping can actually be observed in condensates which 
are realizable at present. 

Finally we want to discuss the problem, to which ex- 
tent the tunneling Hamiltonian is still applicable in a 
regime, where the barrier Vb is only slightly larger than 
the chemical potential, as was assumed above. Now at 
least on the mean-field-level, the quality of the transfer 
Hamiltonian model may be tested by comparing its pre- 
dictions with those obtained from a numerical integration 
of the time dependent GPE in the geometry of Fig. [l]. 
Due to the periodic boundary conditions in the direc- 
tions parallel to the potential barrier, the problem is ef- 
fectively one-dimensional. The time evolution of the con- 
densate wavefunction can thus be easily determined using 
a split-operator Fourier technique. In Fig. ||, the normal- 
ized number of particles which have tunneled through the 
barrier is shown for a case with d — 2£ « 0.3 /jm and a 
comparatively large value A/io = 0.2/2. Time is measured 
in units h/p,, corresponding to about 0.6 msec at t = 30 
for our values above. The height of the potential barrier 
is lowered continuously from an initial value 10/1 
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FIG. 4. Comparison of the normalized particle number dif- 
ference z(t) obtained by a numerical integration of the GPE 
(solid line) and from the semiclassical dynamics based on the 
tunneling Hamiltonian (dashed line) for a high and narrow 
potential barrier with d — 2£ and Vb ~ 3/2 



at t = 10 to Vb = 3/i at t = 20. Charging effects 
are contained intrinsically in the GPE or via the JJv- 
contribution to the chemical potential difference; the ini- 
tial conditions are v(0) — ip(0) = 0. Evidently, the re- 
sult obtained within the transfer Hamiltonian model is 
in excellent agreement with that from a numerical solu- 
tion of the GPE, in the limit of a high and narrow bar- 
rier. Note that there are no adjustable parameters here, 
with all quantities being determined by the microscopic 
parameters of the Bose system. Remarkably, the trans- 
fer Hamiltonian still provides a reasonable approximation 
even for the case of relatively wide and low potential bar- 
riers, which are relevant experimentally. For example, in 
Fig. ||, the function z(t) is plotted for an identical situa- 
tion as above, however with d — 5£ sa 0.8 /an, Vb = 1.4/2 
and AfiQ = 0.1/2 Although the semiclassical dynamics 
based on the transfer Hamiltonian does not reproduce 
the higher harmonics present in the solution of the GPE, 
which are essentially a result of the nonadiabatic lowering 
of the barrier, it still provides a reasonable approxima- 
tion even in this rather extreme case of a barrier which 
is only slightly larger than /2. 
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FIG. 5. Same as Fig. 4, but now for a barrier with d = 5£ 
and Vb = 1.4/2 

VIII. CONCLUSIONS AND OUTLOOK 

In summary, we have developed a microscopic theory 
of Josephson tunneling in weakly interacting BECs. It 
is essentially the analog within Bogoliubov theory of the 
standard work by Ambegaokar et.al. [ p4|j3l| on the tun- 
neling Hamiltonian description of the Josephson effect 
between BCS superconductors. Apart from the fact that 
the Bogoliubov theory only applies far below the con- 
densation temperature T c , while BCS is valid right up 
to T c , there are a number of further crucial differences. 
Most importantly, the Josephson effect in BECs arises 
already in first order in the tunneling amplitude. As a 
result, the Josephson current dominates any other con- 
tributions in the limit of weakly coupled condensates, to 
which our discussion has been restricted. Nevertheless, 
the second order dissipative currents have a strong in- 
fluence on the dynamics of coupled condensates, because 
they remain finite even at zero temperature due to the 
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absence of a gap. Remarkably, the dominant contribution 
to the normal current, which leads to a damping of the 
Josephson oscillations, arises from c-nc tunneling and is 
thus proportional to the condensate density. The explicit 
calculation of the total current has been performed here 
in a very simplified model. It has the advantage, how- 
ever, in providing analytical and physically transparent 
results, allowing to determine the currents in a concrete 
experimental realization from the knowlegde of the bulk 
condensate properties, the barrier height and the effective 
contact area. The fact that tunneling in a Bose system is 
possible at all energies, made it necessary to specify the 
associated matrix elements in much more detail than in 
superconductors. There, only the Fermi energy is rele- 
vant, and thus the matrix elements can be replaced by a 
constant which is fixed by the normal state conductance. 
The effect of 'charging' for condensates with a constant 
total particle number has been shown to be quantita- 
tively important in realistic situations, even though the 
relevant energy U is much smaller than the Josephson 
coupling energy Ej. 

Regarding the prospects for an experimental obser- 
vation of the Josephson effect in condensates of dilute 
atomic gases, our estimates show that this requires bar- 
riers which are only slightly higher than the chemical 
potential. With presently available condensates, the ob- 
servation appears to be possible with 23 Na, not, however, 
with heavier atoms. One of the problems, for instance, 
which may suppress the small oscillating Josephson cur- 
rents, are random fluctuations in the barrier height due to 
fluctuations in the laser intensity, a complication which 
has not been considered so far. 

On the theoretical side, quantitative calculations for 
realistic geometries of coupled traps (in particular be- 
yond the Gross-Pitaevski level) would be useful and also a 
careful consideration of finite temperature effects. For ex- 
ample, the currents due to particles in the thermal cloud 
have been estimated in [fL7[ , however more work needs 
to be done in this direction. Finally, since Josephson os- 
cillation frequencies are typically of the same order than 
those of collective modes [p] of the individual conden- 
sates, it is important to investigate a possible coupling 
between the intra- and inter- well dynamics. 

We gratefully acknowledge helpful discussions with T. 
Esslinger and I. Bloch on the experimental aspects of our 
work. 
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